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Abstract. This paper is concerned with regular flows of incompressible weakly viscoelastic 
fluids which obey a differential constitutive law of Oldroyd type. We study the newtonian limit for 
weakly viscoelastic fluid flows in M N or T for N = 2, 3, when the Weissenberg number (relaxation 
time measuring the elasticity effect in the fluid) tends to zero. More precisely, we prove that the 
velocity field and the extra-stress tensor converge in their existence spaces (we examine the Sobolev- 
H s theory and the Besov-Bj' 1 theory to reach the critical case s = N/2) to the corresponding 
newtonian quantities. This convergence results are established in the case of "ill-prepared"' data. 
We deduce, in the two-dimensional case, a new result concerning the global existence of weakly 
viscoelastic fluids flow. Our approach makes use of essentially two ingredients : the stability of the 
null solution of the viscoelastic fluids flow and the damping effect, on the difference between the 
extra-stress tensor and the tensor of rate of deformation, induced by the constitutive law of the fluid. 

1. Introduction, main results and notations. In this paper we investigate 
the Newtonian limit of weakly viscoelastic fluid flows of Oldroyds'type in = M N or 

n = t n . 

The dynamics of an homogeneous, isothermic and incompressible fluid flows, is 
described by the partial differential derivatives system given by : 

f p(u' + (u.V)uj = / + diver 
[ divu = 

Here p > is the (constant) density and / is the external density body forces, u = 
u(t, x) is the velocity vector field and a — a(t, x) is the symmetric stress-tensor, which 
is split into two parts : a = —pld + r where —pld is the spherical part (p = p(t, x) 
the hydrodynamics pressure) and r is the tangential part or the extra-stress tensor. 
The fluid is called Newtonian if r can be expressed linearly in terms of the rate of 
strain tensor D[u] = i(V« + Vm t ), i.e. 

(1.2) t = 2rjD[u] 

where r\ is the viscosity coefficient of the fluid (in this case (jl.ip is the Navier-Stokes 
system). A fluid for which (|1.2p is not valid is called non-Newtonian or complex fluid. 

Infortunately no universal constitutive law exists for non-Newtonian fluids (see 
for instance [H]). In this paper we consider a class of fluids with memory. For this 
kind of fluids, the extra-stress tensor at a time t depends on D[u] and its history. A 
model taking into account this properties is the Oldroyd's one. The constitutive law 
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of Oldroyd's type [TU] is given by : 

where < A2 < Ai, Ai is the relaxation time and A2 the retardation time. The 
symbol §| denotes an objective (frame indifferent) tensor derivative (see [11]). More 
precisely, 

V T 

-2- = r' + (u • V)t + tW - Wt - a(Dr + tD) 

with W[u] = §(Vit — Vu T ) is the vorticity tensor and a is a real number verifying 
— 1 < a < 1. The limit case Ai > and A2 = corresponds to a purely elastic fluid 
(which is excluded in our analysis), while the limit case Ai = A2 = corresponds to 
a viscous Newtonian fluid. 

The constitutive law (| 1 . 3[) is not under an evolution form. This equation can be 
transformed into a transport equation by spliting the extra-stress tensor into two parts 
t s + t p , where t s corresponds to a Newtonian part (the solvant) and t p to the elastic 
part (the polymer). Setting r s = 2tj(1 — uj )D [u], with uj defined by < uj = 1 — ^ < 1, 
it follows from (|f .3[) that r p satisfies the following transport equation : 

(1.4) Tp + Al ^=2^D[ u ]. 

From now on we shall denote t v — t and rewrite (jl.ip and (|1.4[) by using dimensionless 
variables, we obtain the following partial differential system : 



Re (u' + (u.V)u) - (1 - uj) Au + Vp = / + div r 

d.r,) { divM = infi, 

e(T' + (it.V)T + g(Vu,T)) +r = 2cjD[w] 



where g is a bilinear tensor valued mapping defined by 

g(Vu, r) = rW[u] - W[u]t - o(D[u]r + t D[u]), 

Re = p— and £ = Ai^ are respectively the well-known Reynolds number and the 
Weissenberg number (U and L represent a typical velocity and typical length of the 
flow). It is worth noticing that the Weissenberg number is usually denoted by We. 
Here, since we will make the Weissenberg number tend to zero, we prefer to denote 
it by e. It is crucial to note that when e = 0, (|1.5p reduces to the incompressible 
Navier-Stokes system : 

f Re(v' + (v.V)v)-Av + Vp = f ]nQ 
{ divv = 

On the other hand, from the definition of the retardation parameter we observe that 
u) = 1 — jji/e where < fi < e is given by [i = A2^. Therefore the Newtonian limit of 
(|1 .5|) is actually a limit with two parameters e and /i. To simplify the study we could 
drop a parameter by assuming that the rate fi/e (or equivalently u>) is constant as 
e tends to zero. In this work, instead of doing this, we will only impose an uniform 
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upper bound on to (= 1 — /i/e) with respect to e. 

System (|1.5[) is completed by the following initial conditions 



(1.7) 



«|t=o = ""o and T| f=0 = t 



Our approach is quite general and uses the two following ingredients : 

• The stability of the null solution of (jl.5|) for a fixed e (see 2, on JR^ or T N and [7] , 
[5], [9] for the case of a bounded domain). 

• The damping of factor 1/e on the quantity r — 2cjD[u] induced by equation (|1.5[) q. 
Our results in the Sobolev spaces are valid for Q = IR N or T N but to simplify the ex- 
pository we will only consider O = IR N and give the necessary modification to handle 
the periodic case. 

The main idea is to cut u and r in low and high frequencies at a level depending on 
1/e. Roughly speaking, forgotting the nonlinear terms, the high frequency part of 
u — v (v is the Newtonian solution, see (|1.6p associated with the initial data uq) will 
satisfy the homogeneous system linearized around the null solution plus a non homo- 
geneous part containing a high frequency term of v. But by the Lebesgue monotone 
convergence theorem, this term will tend to zero in the appropriate norms. The sta- 
bility of the null solution (cf. [2], [7]) will then force the high frequencies of u — v and 
e 1 / 2 r to remain small (recall that (u — v)(0) = 0). On the other hand, the remaining 
frequencies will tend to zero due to the damping effect on r — 2wD[m] which we will 
use in the same time as a smoothing effect. We will describe the main steps of the 
proof in Section 1.3. 

Note that our analysis is in the spirit of numerous works on the incompressible limit 
of compressible Navier-Stokes equations (see for instance [3] and references therein). 
However, our analysis is in some aspects easier since there is a damping effect relating 
to the small parameter whereas in the incompressible limit it is a dispersive effect. 

In our knowledge, no such result exists in the literature concerning our study, 
i.e. the newtonian limit of non-newtonian fluid flows. Moreover, our global existence 
result for regular weakly viscoelastic fluids flow in dimension 2 (see Corollary II. ljl is 
new and, in particular, not contained in the global existence results of [2]. 

1.1. Function spaces and notations.. In the sequel C denotes a positive 
constant which may differ at each appearance. When writing x ^ y (for x and y 
two non negative real numbers), we mean that there exist C\ and Ci two positive 
constants (which do not depend of x and y) such that C\x < y < C2X. When writing 
x < y (for x and y two non negative real numbers), we mean that there exists C\ a 
positive constant (which does not depend of x and y) such that x < C\y. 
V will denote the Leray projector on solenoidal vector fields. 
For 1 < p, q < 00, we denote by || • \\lp the usual Lebesgue norm on Q = M N , 




and by || • the space-time Lebesgue norm on ]0,t[xf2, 
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with the obvious modification for p,q = oo. For s G Ft, we denote by || ■ the usual 
Sobolev norms on £1 = 1R N , 

\Mh. = ([ (i + \er\m\ 2 ^) 1/2 

where v is the Fourier transform of v. The corresponding scalar product will be 
denoted by ((•, -))h s ■ Finally, for any e > we introduce the following Fourier 
projectors 

(1.8) ^(0=X[o )E «](l£l)/(0 and QJ($ = X] e -,oo[(lfl)/(0 , 
where a > will be specified later. 

1.1.1. Homogeneous Besov spaces. Let ip in S(IR) such that $ is supported 
by the set {z / 2 _1 < |z < 2 } and such that 

(1.9) 5^(2-'ar) = 1 , * ^ . 

ieZ 

Define <p by 

(1.10) £ = 1-^(2"^) , 

and observe that y> G T>(JR), if is supported by the ball {z / \z\ < 2 } and (p = 1 
for |z| < 1. We denote now by Aj and Sj the convolution operators on 1R N whose 
symbols are respectively given by $(2 — and <£(2 _J '|£|) where £ G JR^ and |£| = 
+ ■■ + £,n- Also we define the operator Aj by 

Aj = A rf _i + Aj + A i+1 , 

which satisfies, 

Aj o A, = A, . 

For s in 1R, the homogenous Besov space B^ilR 1 ^) (to simplify the notation we 
will simply denoted it by B S (IR N ) ) is the completion of S(]R N ) with respect to the 
semi-norm 

(1-H) = ||{2^||A j (/)|| L2 }|| il(Z) . 

1.2. Main results. Theorem 1.1. Let N = 2,3 and let (uq,t ) G H S (R N ) x 
H S (]R N ~) and f G L 2 ^ oc (M; H 3 ^ 1 ) with s > N/2. Let v be the Newtonian solution 
satisfying hi. 6}) with initial data uq and let < To < oo such that v G C([0, To]; H s ). 
Then for any S G]0, 1[ there exists 

eo = e (N, Re, 5, \\v\\ l ~h; \\Vv\\ l1qH s, \\t \\ h ., \\Pf\\ L ^ H .-i) > 

such that for any < e < Sq the system \1.5\) . with 



(1.12) 



< oj < 1 - S, 
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admits a unique solution 

u £ eC([0,T o ];H s ), Vu £ eL 2 (0,T Q ;H s ), r £ e C([0,T }; H s ) . 

Moreover, 



(1.13) u e — > v inC([0,T ];H s ), 

(1.14) t £ - 2ujT>[u £ ] — > in L 2 (0, T ; H s ) , 

(1.15) e 1/2 r £ — ► in C([0,T o }; H s ) . 



Recalling that in dimension two, the solution of the Newtonian problem exists for 
all positive times, we deduce the following result. 
COROLLARY 1.1. In dimension 2 there exists 

e = £ ( Re, S, \\v\\ LSH ; ||Vv|| ZSo ff., WMh; H^/IU^-O > 

such that for any < e < £q the solution of \1.5\) given by Theorem \1.1\ exists for all 
positive times. 

Remark 1.1. Note that Theorem \1.1\ is a convergence result for "ill-prepared" 
data. Indeed the quantity tq — 2luJD[uo] is not assumed to be small with e. Moreover, 
this is a singular limit result since r and D[u] do not belong to the same function 
space. In particular, D[uq] is not as the same level of Sobolev regularity as tq. 

Remark 1.2. According to the introduction, the Newtonian limit process is ac- 
tually a limit process with two parameters e and fi tending to zero with < fi < e. 
The assumption (l.lfy) of Theorem \l.l\ means that we impose the following additional 
conditions on the rate life as (e, /i) tends to zero (0,0): 

<5<^ = ^<1 
e X\ 

for some fixed 1 > 5 > 0. 

As mentioned in the introduction, the use of Besov spaces enables us to reach the 
critical index s = N/2. 

Theorem 1.2. Let N — 2,3 and let (u ,t ) £ B N / 2 ~ 1 (M N ) x B N ^(M N ' 2 ) and 
f G Lj oc (_B JV / 2_1 ). Let v be the Newtonian solution satisfying il.6\) with initial data 

Uo and let < To < oo such that v s C([0, To]; B N ^ 2 ~ X ). There exist < ujo < 1 and 
£o = £o(N, Re,uo, ||to||bjv/2,'P/,Uo) > such that for any < s < Eq the system 
il.5\) . with < oj < ujq, admits a unique solution 

u £ e C{[0,T Q ]:B N / 2 - 1 ), u £ G L\0,T a ;B N ' 2+1 ), r £ G C([0, T ]; B N ' 2 ) . 

Moreover, 

(1.16) u £ — > v m C7([0,T ];B Ar / 2 - 1 ) , 

(1.17) t £ ~2luT>[u £ ] — > in L\0,T a ;B N/2 ) , 
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(1.18) £ 1/2 r £ — » mC([0,r ];B N / 2 ). 

In dimension two, using the classical global existence result in B° for the Newto- 
nian problem (see for instance [3]), we get the following corollary. 
COROLLARY 1.2. In dimension 2 there exists 

e = e ( Re,uJ , \Hl^b», W^Wl^, \\to\\bU WPfh^B") 

such that for any < e < Eq the solution of M.5\) given by Theorem \1.2\ exists for all 
positive times. In the surcritical case, s > N/2, we get similar results by considering 
non homogeneous Besov spaces. Note that £o depends then explicitly on some norms 
of v. 

Theorem 1.3. For s > N/2, Theorem [T7B and Corollary ] also hold by replac- 
ing the function spaces B™/ 2 " 1 by B^^B^ 2 " 1 and B N / 2 by B s n B N / 2 . Moreover, 
Sq will depend now explicitly on some norms ofv andVf. More precisely, for s > N/2, 
we have 

£o = e (N, Re,u , \\v\\ l ^ b n/2-i,\\Vv\\ l ^ B s,\\t \\ b n/2, \\Vf\\ L i b n/2-i) . 



1.3. Sketch of the proof of Theorem 11.11 In this subsection we want to 
explain the main steps of the proof of Theorem 11.11 Note that Theorems 11.21 follows 
from the same arguments. To simplify we drop the nonlinear terms in (jl.5p . The first 
step consists in noticing that W := u — v satisfies the following system: 



(1.19) 



Re W t - (1 - uj)Q e AW - P e AW = P e P(divr - luAu) 
-ujQ e Av + Q e VdbfT 

div W = 

er t + Q e r = 2ujQ e T>[W] + 2uQ £ T>[v} - P £ {t - 2wD[u]) 



where P £ and Q e are the projectors on respectively the low and the high frequencies 
defined in fLS]) . 

Projecting on the high frequencies with Q e (see (|1.8|) for the definition), proceeding 
as in [2 j , it is easy to check that we get a differential inequality close to 

j t (\\Q s W\\ 2 H s+e\\Q 5 rf m ) + \\Q B VWf H . + \\Q s r\\ 2 H , < \\Q e Vv\\ 2 H s 

where we drop all the constant to clarify the presentation. Therefore, since W(0) = 0, 
e — > and, by the Lebesgue monotone convergence theorem, ||Q e Vi;|| L 2 Hs — > 0, we 

infer that ||Q e W^||i5? H* goes to zero with e. Now, to treat the low frequency part, 

we observe that, computing P g f (|1.19|h j^-D[ (|1.19|) i] and taking the i? s_1 -scalar 

product of the resulting equation with Z := r — 2wD[u], we obtain something like 

(1.20) Illi^U^-* + \\\P e Z\\ H ^ < \\P e r\\%. + \\PJWls-, . 

On the other hand, P„ (fTT9)h can be rewritten as 

eP e T t + e p P £ r = 2we /3 P e D[Fy] + 2loe i3 P £ T>[v] - (1 - e fi )P e Z , 
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where < /3 < 1 will be specified later. Therefore, taking the £P-scalar product of 
this last equality with r and adding with the scalar product of (|1.19[) i with W we get 
a differential inequality close to 

Adding this last inequality and e 2 ^ ()1.20p we finally obtain 

J t (\\PeW\\%,+s\\P £ rf m ) + \\PeVWf Hs +e%P e rf H . + ^\\P e Zf H ..r 
<£nPjf H .-i+e?\\P e Vvtf H . 

since e~ p \\ P s Z\\ 2 Hs < e^lP^Zf^^ < ^^H^H^.-i as soon as 1 - 3/3 - 2a > 0. 
This last inequality enables us to conclude for the low frequency part. Note that we 
used the damping effect also as a smoothing effect. 

2. Proof of Theorem [TTTJ. Let us recall the following existence theorem proven 
by J.-Y. Chemin and N. Masmoudi [2]. 

Theorem 2.1. Let (u ,t ) e H S (M N ) x H s (M n2 ) with s > N/2. Then there 
exists a unique positive maximal time T* and a unique solution 

(«,r) S C{[0,T*[;H s )nLf oc {0,T*;H s+1 ) x C([0, T*[; H s ) 

Moreover, if T* < oo then V N/2 < s' < s 

(2.1) limsupflKt)!!^ + \\r(t)\\ HS A = +x> 

Remark 2.1. Actually in '21 the following sharper blow up condition is derived 

T*<OQ=> [ \\Vu(t)\\ La * + \\T(t)f L ao dt = +<x> , 

Jo 

but for our purpose the classical blow-up condition will be sufficient. 

Let us also recall a commutator estimate and classical Leibniz rules for fractional 
derivatives. 

Lemma 2.2. Let A be the Laplace operator on M , N > 1. Denote by J s the 
operator (1 — A) s / 2 . 

• For every s > N/2, 

(2-2) \\[J s J]9\\L*(n») < l|V/||H. (K * ) ||ff|| H .-i (Jl *) . 

• For every s > 0, 1 < q, q' < oo and 1 < r, p,p' < oo with l/p+ 1/q = l/p' + l/q' = 
l/r, 

(2.3) || J S {fg)\\ L r(lRK) < \\J S f\\hP(]R N )\\9\\Li(lR N ) + \\f\\Lv'(Fi N )\\J S 9\\Lp'(m N ) ■ 

• For every p, r, t such that r,p > t, r + p > and r + p — t > N/2, 

(2-4) Wfg\\Ht(lR N ) S Wf\\HP(lR N )\\g\\Hr(lR") ■ 
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Proof. (|2.3p and (|2.4|) are classical and can be found in [5] and 0. (|2.2p is a 
variant of Kato-Ponce's commutator estimates. It is proven in [13] in dimension 1 but 
the proof works also in dimension 2 and 3. 

To treat some nonlinear terms in dimension 2 we will need moreover the following 
Gagliardo-Nirenberg type inequality (see for instance [1]). 

Lemma 2.3. Let N > 2, for u 6 H 1 (1R N ) the following Sobolev type inequality 
hold for any 2 < p < +oo such that i — < ~ : 

( 2 -5) II"IUp(JR w ) ^ IMIl^R") H Vu ll ■ 



2.1. Estimate on W = u — v and e 1 ' 2 t.. We start by deriving a differential 
inequality for the ff s -norms of W and £ 1 ^ 2 r. The high frequency part of this inequality 
is directly inspired by the stability proof of the null solution in [2]. This will enable 
us to control the very high frequency part (Q £ u, Q £ t) of the solution. The other part 
(P e u,P e r) will be treated by using the damping effect. 

For e > fixed, Theorem 12.11 gives the existence and uniqueness of the solution 
(u e , t £ ) of (TO]) in C([0, T*[; H s ) n L 2 oo (0, T*\ H s+1 ) x C*([0, T*[; H s ) for some T £ * > 0. 
To simplify the notations, we drop the index eona and r in the sequel. Setting 

Z = t- 2wD[u] 

we have the following estimates : 

Lemma 2.4. For e > small enough, the solution (u,r) of hi .5}) satisfies for all 
<t<T E * and < (3 < 1, 

+ £lMI^) + \\\P S VW\\% S + ^l\\ Qs vw\\% s 



4e"' 3 



< (1 + Jll^llff. + MQeVv\\%. + 8<J£ \\P e Vv\\ 2 H3 



( 2 - 6 ) +7?-^(l|V«[||. + IIVHI^OII^II^ + ^^IIVull^Hlr. 

(1 — w 

whenever < < 1. Moreover, for < uj < 10~ 2 , if ZioZrfs 

d ( Re m*ru2 , e„ ||2 \ , 1„ D „ W „||2 , (1 - w )„^ V7WII2 



+JIIQ £ * + ^WI&. 

< (l + 4e- /3 )||P £ Z|| 2 fs + 8w 2 ||Q s VHIh« +8u:' 2 e f3 \\P e Vv\\ 2 H . 

(2.7) + ^^(||V U || 2 ffa + ||V^| 2 f 0ll^ll^+C e 2 ^||V U || 2 Hs ||r|| 2 , s . 

(1 - uj) 2 

Proof. Notice that W verifies the equation 

Re (W t + P(u.W)w) - AW = Pdivr - luAu - Re V{W.V)v 

= P s "P(divT - wAuJ - uQ £ Av 

(2.8) +Q £ Pdivr - ujQ £ AW - Re P(VK.V)v . 
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Therefore, multiplying scalarly flXS]) by W in H S (M N ), using Cauchy-Schwarz, Lemma 
and that u is divergence free, we obtain 

\ Re±\\W\\ 2 Ha + HPeWH^ + (1 - u>)\\Q e VW\\ 2 Hs 

< {(Q £ dWT,W)) H s + \\P e Z\\ H 4VW\\„s+ij\\Q E Vv\\ H 4Q e VW\\H* 



(2.9) 



+C Re 



((J s (u.V)W, J S W)) L 2 



+ Re 



((J'(W.V)v, J S W)) L 2 



To estimate the second to the last term of (|2.9j) . we rewrite it with the help of a 
commutator and apply Cauchy-Schwarz to the term containing this commutator to 
get 

((J s {u.V)W, J S W)) L 2 < (((u.V)J s W, J S W)) + [J s ,{u.\7)]W \\J S W\\ L 2 . 

L 2 

Since u is divergence free, the first term of the right-hand side of this last inequality 
cancels by integration by parts. Estimating the second term thanks to Lemma 12. 2\ 
we then obtain 



((j s (u.v)w,rw)) L 2 



< CHVul 



\vw\ 



\w\ 



Now, to estimate the last term of the right-hand side (|2.9p we have to distinguish the 
cases N = 2 and N = 3. 

• N = 3. Then by Lemma [2721 Holder, Sobolev and Young inequalities, we get 



((J s (W.V)v, J S W)) L 2 < \\J s {W.V)v\\ L6 /4J s W\\ L6 



< \\J s w\ 



L 2 



|V«| 



\w\ 



J s Vv\\ T ,2 \\\J S VW\\ T ,2 



<\\W\\ Hs \\Vv\\ H s\\VW\\ H s ■ 
N = 2. In this case, using Holder and Lemmas 12.2112.31 we infer that 

((J s (W.V)v), J S W)) L 2\ < ||J'(W.V)i;|| £ s/ 3 ||J'W|Ua 

< (||J S H/|| L6 ||V«|U2 + \\W\\ L e\\J s Vv\\ L 2^\\J s W\\ La 

< \\J s W\\ L 4J s Vv\\ L 2\\J s W\\ L3 

< II Wllja-- 11 Vwlljy. II VWllir- • 



By Young inequalities it thus follows from (|2.9p that 



Re d 
~~2~~dt 



\Wf Ha + -\\P e VW\\ 2 H , 



(1- W ) 



< ((Q £ divr, W)) H . + \\P e Z\\ 2 H . 



4(1- u) 



\\Q>eVW\\ 2 Hs 
\QeVv\\ 2 HS 



(2.10) 



C Re 



(||V u || 2 ffs + \\Vv\\ 2 Hs )\\W\\ 2 Ha 



On the other hand, for < (3 < 1, observing that 

t - 2wD[u] = Q £ t - 2ujQ £ (T>[W] +T>[v}) 

+ (1 - e f) )P £ Z + e f3 (p e T - 2ljP £ {T>[W] +~D[v]) 
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we deduce from (|1.5|) that r satisfies the equation 

e(r t + (u.V)t + g(Vu, r)) + Q £ t + e f3 P £ r = 2ujQ e T>[W] + 2wQ e D[v] 
+2uie i3 P e D\W] + 2uje 13 P E T)[v] - (1 - e p )P e Z . 



Taking the H s scalar product of this equation with t, using Lemma 12.21 Cauchy- 
Schwarz and Young inequalities we get 

ll l|r||2ffs + \^ T ^ + y 11^* ^ M(QMW],t)) h . +8e-P\\P £ Z\\l., 
+8uj 2 \\Q e Vv\\ 2 Hs +8u 2 e f, (\\P e VW\\ 2 H . + \\P e Vv\\ 2 H .) 

(2.11) +Ce 2 -e\\Vu\\ 2 Hs \\T\\ 2 Hs . 

We now separate the two cases : 

• uj ^ 0. Then, adding (|2 . 10[) and (|2. 1 ip / 2a; we notice that the first term in the 
right-hand side of (|2 . 10[) and (|2 . cancel each other and (|2.6[) follows. This gives 
(|2.6|) for e small enough since /3 > 0. 

• < uj < 10~ 2 . Then adding (f2~T0]l /2+ ((2TTT]) . estimating the two remaining H s - 
scalar products by integration by parts, Cauchy-Schwarz inequality and Young in- 
equality, one obtains (|2.7[) 

2.2. Estimate on Z = t — 2wD[u]. We will now take advantage of the damping 
effect on Z = t — 2cjD[u]. 

Lemma 2.5. The solution (u,t) of \1.5\) satisfies for all e small enough and 
0<t< T*, 

(2.12) +t^-( « e ll v »llH- + IMIh.)IMI«. ■ 

1 — UJ 

Proof. We apply -p^-D[-] to (|1.5P i and substract the resulting equation from 
(11.5^ 9 to obtain 



(2.13) Z t - ^r^AZ +-Z = -A - f 2 

Re e 



where 



and 



h = ^D[Pdivr] - ^^Ar + ^0[Vf\ - 2uT>[T{u.V)u] 



f 2 = V(u.V)t + g(Wu, t) . 
Taking the _ff s ~ 1 -scalar product of (]2 . 13[) with Z we get 

Re (1 — uj) Re (1 — uj) 

(2-14) +^i||( u .v)u||^_ l + 4||(«.V)r||^_ l + ||g(Vu, r)f H ^ 
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where we used that 
2w 
~Re 



2uj 

{(r>[PdiVT],Z)) H s~t < C—\\diVT\\ H .-x\\VZ\\ B .-x 

<—\\vz\\ H ^ + c- 2 



Re " ' - ' Re(l-ct)) 

Finally to control the nonlinear terms we notice that thanks to (|2.4[) , 

Ho.V6HiT.-i < llollH.llV6Hir.-i 
which concludes the proof of (|2.12[) . 

2.3. Convergence to the Newtonian flow. We give here the proof in the 
case 10~ 2 < bj < 1 — 5. The case < u> < 1CP 2 is simpler and can be handled in the 
same way by using (|2.7p instead of (|2.6[) . 
Adding (EU) and e 213 jUg) . we obtain for e small enough 



-28 



K^Oir. + ^IMIi. + Vll^ 

+ (±=A\\ VW f Hs + ±-(\\Q e r\\l,+e^\\P e r\\l, 



.2/3-1 



-\\Z\\ H s 



< 8w 2 ||Q £ Vw||^ a + 8u)' z e p \\P e Vv\$i. + Ce 1_/J ||Vit[[ H 



4w 



+ 7?^(I|V«||^ + ||V W ||L 



(l-u;)* 



(2.15) 



-C- 



z 2 < 3 r 1 



Re 



(1- W ) 

Here, we used that for e small enough 

(1+^) 2 
Re(l-w)' 



||P/|| 2 ff3 _ 1 + (Re||V U || 2 ffs 



-||^.)II«IIh. 



2/3 | 



Ik 



8w 



Ik 



and 



(1 



\P s Zf H s < {l + 2Q 2 e-P)e- 2a \\P £ Z\\ 2 Hs ^ < 



.2/3-1 



\\Z\\ H s 



as soon as /3 > and 1 — 3/3 — 2a > 0. 

From now on we thus take to simplify (a,f3) = (1/8, 1/8). Setting 



Re f 
X.(t) = —\\W(t)\\%. + —\\r(t)\\ 3 H . 



r 2/3 



-\\Z\\ H s-l 



(1-0,) 



1 



e 



4 l|VW|| H3 + -||Q £ r|| Hs + -||P £ r|| Hs , 4 

we infer that X s satisfies the following differential inequality 

d e 2(3 

-X. < 8u 2 \\Q E Vv\\ 2 Hs +8co 2 s (3 \\P E Vv\\ 2 Hs +C-— ; 

at Re (1 — uj 



2/3-1 



■\\Z\\h-*d8 



■WPfWl* 



C(Re,5) e 20 \\tW hb + \\VuW H s + \\Vv\ 



x s 



(2.16) +C 



(1- W ) 



(Re^||V U || Hs +^||r|| Hs )ll«llH. 
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where we rewrite u as W + v and use the triangle inequality when necessary. Hence, 
Gronwall inequality leads to 

X s (t) <exp[c( Re,£)(>||T||| ?Ks + ||V«||£ ?HJ + IIV^ 2 ^ 

-2/3 



X,(0) + 8uj 2 \\Q £ \7v\\% Hs + 8e^||V«||| fff . + C—\\Vjf m .-, 



(2.17) 



— ( Ree^llVi 



\r\\% Ht )\\v 



Rewriting u as v + W, we finally obtain 

<exp[(7( Re,<5)(||V«|| 2 ?HS +Jf s (t) 



X s (0) + 8^ 2 ||Q E V W ||^+8 e ' 3 ||V U ||^ 3 +C'^ i7 ||P/||| |HS _ 1 



(2.18) 
where 
(2.19) 



C 



(1+ Re)e /3 



W||ioc ff .X.(t)+C 



Ree 2 ' 3 



iVul 



e e 2 * 3 
X s (0) - ^llroll^ + — ||r - 2 W D| Uo ]||^- 2 



Let us now assume that T* < Tq. Since (|2. 18|) holds for any N/2 < s' < s, noticing 
that 

\\QeVv\\ L 2 H ., <S a ^\\Vv\\ L2tHs , 

we deduce from (|2.18p . (|2.19[) and the continuity of t ^ X s /(t) that there exists 
£o(s, ||7o||i?», II«o||h=, ||Vw|| L 2 Hs , \\v\\l^ H ') > such that for any < e < e and 
any < t < T* 

(2.20) X s >(t) < C £ ™»(/3,2 Q ( S -s')) < Ce min(l,s- S ')/S 



which contradicts (|2.1j) of Theorem 12.11 This ensures that T* > To. Now, since by 
Lebesgue monotone convergence theorem 



(2.21) 



\QeVv\\ L 2 H s— >0 



it follows from (|2~T8l) - (|TT9| that X S (T ) -> as e -> 0. This proves (fL~13|) and 
(|1.15|) . To prove (|1.14|) we observe that from this last limit and (|2.2ip . ||Q e (r — 
2cjD[u])|| L 2 ffs -> and e 2 ' 3 - 1 ||P e (r - 2o;D[u])|| 2 A -► 0. This yields the result 



by Bernstein inequality since 2(3 — 1 + 2a < 0. 



2.4. The periodic setting. Let us give here the modifications needed to handle 
with the case n = T N , N = 2,3. It is worth noticing that Lemma [2~2l holds also with 
fl = T N . On the other hand, the Sobolev inequality (|2.3p does not hold for general 
functions in T N but holds, for instance, for zero-mean value functions. Note that if 
f(t) has mean value zero for all time t > then using the invariance by Galilean 
transformations, u > u(t,x — zt) + z with z £ M 3 , we can assume that u has zero 
mean- value for all time and we are done. Otherwise, we have only to care about the 
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treatment of the nonlinear term (W.V)v in (|2.8p . Denoting by W_ the L 2 -projection 
of W on zero mean-value functions, we rewrite (W.X7)v as 



(2.22) 



(W.V)v = (W.V)v + ( ( W)Vv 
Jo. 



We take the .fP-scalar product of (|2.8|) with W_ and add with the L 2 -scalar product 
of (12. 8j) with The if s -scalar product coming from the first term of the right-hand 
side of (|2.22j) can be treated as in M N . For the second term, we observe that 



w\ 



((Vv,W)) H , 



< \\W\\ l 4Vv\\ h .\\VW\\h-i 
On the other hand, concerning the L 2 -scalar product we notice that 



(((w.v)«, 



w 



L 2 



[[{W.V)v,W^j l2 + (J W) ([Vv,W 



L 2 



<\\w\\ L s\\vv\\ L 2\\w\\ L3 + \ / w\\\Vv\\ l 4w\\l* 

< \\W\\ H s\\VW\\ L 2\\Vv\\ L 2 

We thus obtain exactly as in (|2.10p . 



Re d 
~~dt 



(llffll^ 



\W\ 



L 2 



< {{Qsdivr,W)) H s + \\P e Z\\% 
C Re 



|V«||^ 



4(1 - w) 



Q> £ VW|lff. 

IQeVi;|| 2 ffs 



The remainder of the analysis is now exactly the same as in M 

3. Proof of Theorems 11.21 and 11.31 In this section we prove a convergence 
result in the Besov spaces B%~ ' , s > iV/2. It will require a smallness assumption on 
the retardation parameter u> but on the other hand will enable us to reach the critical 
regularity space for (| 1 . 5j) . Note that our smallness assumption on the retardation 
parameter is the same as the one in [5] to get the stability of the null solution in such 
function spaces. 

Let us recall the following well-posedness result derived in [2] . 

Theorem 3.1. Let(u ,T ) e B s - 1 (M N )r\B N ^- 1 {M N )xB s (M N2 )nB N / 2 (M N2 ) 
with s > N/2. Then there exists a unique positive maximal time T* and a unique 
solution 

(«, t) e C([0, T*[; B 3 - 1 DB N / 2 - 1 )nLl c (0, T*; B s+1 nB w/2+1 ) x C([0, T*[; B S C\B N ' 2 ) 
Moreover, if T* < oo then 



(3.1) 



limsupf ||u(£)|| B iv/2-i + \\r(t)\\ B N/2 

t/T* v 



-oo 
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We will make use ol the following classical commutator and product estimates 
(see for instance [2], [3] and [H]) 

Lemma 3.2. For all s e]l - iV/2, 1 + N/2[ we have 



(3.2) ||A J -[(o.V),A j ]6|| a < 2-^- 1 ) 7i ||Vo|| BJ » /a+ i||6|| B .-i, 
mth ll7illii(Z) ~ L 

i*br aZZ si, S2 < N/2 with s\ + S2 > if /loZcfe 

(3.3) \\ab\\ Bsi+S2 - N/2 < ||a[[ B » 1 ||6[[ B » 2 . 



For any e > we divide Z into the three following subsets 
I :=%*_ = {j e Z, < V < 1}, J e := {j e Z, 1 < V < eT"} and JT £ := {j e Z, 2 J > e""} 
and for any subset N C Z we denote by || • the semi- norm 

3.1. Estimate on W and er. Lemma 3.3. The solution (u,et) of il.5)) satis- 
fies for allO <t<T* 

j t (Re\\W\\ B ,-i+4e\\T\\ B s) 

+ [(1 - w)/2 - 16^11^11^+! + ||W|| B| +i £ + 2||r||^ + 2eP\\T\\ B}uJe 
< 5\\Z\\ Bfujc + IMMIbjH + 16^(|| W\\ B ^ + IHb^J 

(3.ty-C e(JLx\\u\\ B N/3+l\\T\\B> + C(\\U\\ B N/2+l + \\v\\ B N/2+l)\\W\\ B s-l . 

Proof. Applying Aj to (|2.8p we have for j e J e , 

Re (d t AjW + V{u.V)AjW^ - (1 - w)Aj-AW = -uAjAv + AjPdivr 

(3.5) + ReA j V[(u.V),A j ]W + Re AjV(W.V)v 
and for j £ I, 

Re (fltAjW + P(u.W)AjW) - AjAW = A 3 Z 

(3.6) +ReA j V[(u.V),A j ]W + ReA/P(W.V)u. 

Taking the scalar product in L 2 (M N ) ol (|3.5p with AjW, using that W is divergence 
free and Cauchy-Schwarz inequality a we get 

i Re^llAjWIll + a-wJIlVAjWIll < ||A,^|| 2 (u>||A 3 At>|| 2 + ||AjdivT|| 2 

(3.7) + Re ||A i P[(«.V), &j]W\\ 2 + Re HAjPpV.V)*)^) . 

We use now that, according to Bernstein inequality ||VAjW||2 > 2- 7-1 1| AjW[[2 and 
divide (|3.7|) by |[AjW[[2. Then, estimating the commutator term thanks to (|3.2p 
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and the last term thanks to (|3.3|) with Si = s — 1 and s 2 = N/2, using Bernstein 
inequalities, it follows that 

Re^\\A J W\\ 2 + ^-^2^\\A 3 W\\ 2 < 2||AA,t;|| L2 + 2|| A,divr|| 2 



(3.8) + 7j 2 J (s 1) (ll"lli3«/2+i + ||v||i3«/2+i)||^||B S -i 

with ||(7j)||ji(Z) ~ I* Multiplying by 2^ s_1 - ) and summing in j e K e , it follows that 



(3.9) < (N| B */2 +1 + \\v\\ BN/ ^)\\W\\ B s-x . 



^-\\W\\ m -i + ±-—L\\W\\ Bl r - 2\\t\\ b . 



Proceeding in the same way with (|3.6p but summing in j S / U J e , we obtain 
(3.10) < (llullB^+i + llwllB^+OHWllfl.-!. 



Now, for j S Z, we infer from (| 1 . 5[) that 

ed t AjT + e(u.V)AjT + Ajt = 2wAjD[u] 

(3.11) -e[(«.V), A,-]t + eA,g(Vu, r) . 

Rewriting A 3 (r - 2wD[u]) as Aj(t - 2uiT>[W] - 2ujB[v]) for j e if e and as 

Ajt - 2u;e' 3 A j {T>[W] + D[«]) + (1 - e ,3 )A J Z 
for j eJU J £ , similar considerations as above lead to the two following inequalities 

(3.12) +4w||u|| s a +1 + Ce \\u\\ b n, 2 +i\\t\\b° , 



and 



(3.13) +4w^|M| B a +1 +Ce||«|| B » /2+1 ||r|| Ss 

-TUJ e 

Adding (0 + (I3~TUD + 4f (l3"T5D + (l3"T2l h dOl l follows. 
3.2. Estimate on r - 2o>D[it]. Lemma 3.4. 



-||Z|| Br2 + -\\Z\\ B .-, < H^||t|| b . + \\Pf\\ B . 



1 \\7\\ <( 1+w ) 

(3.14) +C* a ln( £ - 1 )(|| U || BN/2+1 + ||r[| B * /B ) |]u|| 



^ll^ + i|mB f <^IH| B . + M | JI .. 

(3.15) +c(\\u\\ b n/2+i + ||r|| B «/ 2 ) . 
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Proof. Applying Aj to (|2.13p and taking the L 2 -scalar product with AjZ we get 

|||A^|| i2 + il^2^||A^|| i2 + i||A^|| i2 < IIA.MU. + \\*jh\\v , 
where 

2uj (1-u) 2uo 

h = — Dpdivr] - [ —^ 1 &t + - jfc D[7>/] - 2ujT>[V(u.V)u] 

and 

f 2 = V{u.W)t + g(Wu, t) . 
Multiplying this inequality by 2^ s ~ 2 \ summing in j S J £ we infer that 

(3.16) +||(u.V)«|| B ,-, + ||(«.V)t||„-. + ||g(V«,T)|| B .-> . 



For s > 1 we estimate the nonlinear term thanks to (|3.3[) with respectively (si, S2) = 
(s - l,iV/2), (s - I, N/2 - 1) and (s - 2, N/2). For s = 1 ( of course N = 2) we 
estimate the first nonlinear term in the same way and use the following lemma to 
estimate the two last ones. This lemma follows directly from the definitions of / and 
J £ and the fact that, for \s\ < N/2, the usual product maps continuouslj0 B^ 8 ' 1 x B s ^ 
into B~ N / 2 '°° (see for instance [H]). Note, in particular, that \J e \ < aln(e _1 ). 
Lemma 3.5. For all si, S2 < N/2 with Si + s 2 = it holds 

(3.17) \\ab\\ R - N/ 2 <a\n(e- 1 )\\a\\ Bsl \\b\\ B s 2 . 
and 

(3.18) \\ab\\ B - N/ 2 + 2 < ||o||fl.i||6|| fl . a . 

We apply this lemma with (si, S2) = (0, 0) and (—1, 1) for respectively the second 
and the third nonlinear term of (|3.16[) to complete the proof of (|3.14p . Finally (|3 . 1 5|) 
can be easily obtained in the same way by using that ||a||s s < ll a llss' for s' < s and 

(EH. 

3.3. Convergence to the Newtonian flow. From now on we set j(u>) — 
(1 — uj)/2 — 16w and assume that < to < u)q with 7(0^0) > 0. 

We proceed as in Section 2.3. For < (3 < 1, we add ^ and g^f fST!)) + (j3~T5|) ) 

to get 

j t ( Re \\W\\ B .-i + 4e\\t\\ B ' + ^(\\Z\\ BZ , + \\Z\\ Bf 

=.2(9-1 , . 

+7 ^o)\\W\\ B s +1 + \\W\\ B ^ +2\\r\\ Bk +^l|r||B fuJ + -5 -(ll^L— + \\Z\\ B .) 
< 16w|Mlu£* + 16we^||«|| B .+i < +e a/J ||^/|| B .-i e 

+Cs||u|| b w/3 +1 ||t||b' +C(||lt|| B JV/2+i + ||u|| B w/a +l )||W|| B .-l 

(^£^ln( £ - 1 )(||t t || BW/2+1 + ||r|| BN/2 ) ||«|| fl ._ 1 . 



^Forl<p<oo, \\f\\ B s, P = ||{2^||A,(/)|| i2 }|| !p(Z) . 
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Here we used that for e small enough, e 13 < min(167(wo), -^r-), £ 2l3 ~ 1 /2 > 5 and 



■>\\z u, <e-^\\z\\ B y l >< e —;z:. l 



4 

2/3-1 



as soon as 



(3.20) < a < 1/2 and < 2/3 < 1 - 2a . 

From now on we set (a, (3) — (1/8, 1/8) so that (|3.20|) is satisfied. Setting 
X s (t) = Re||W(t)|| B .-i +4e||r|| Bs + e 2p \\Z\\ B .-2 

IUJ £ 

+ f ^\\W\\ B . +1 + (||r]]B ie +^]]r||B» UJe ) + ^-(\\Z\\b--> + \\Z\\B f ) da, 
we infer that 

j t X s {t) < 16cj\\v\\ b s kV + l&^IMI^ +e 2f3 \\Vf\\ B s- k 

+C( Re,uj) [\\W\\ BN/2+ i + \\v\\ BN/2+ i +e 2f) \n(e- 1 )\\T\\ BN/2 ~]x s 
+C*ae 2/3 ln(£- 1 )(||t)|| B »/2+i + \\W\\ BN/ 2 + i + \\t\\ b n /3 ) \\v\\ b ,-i . 
By Gronwall lemma we infer that 

X s {t) < exp(c(u, Re)(||«|| i i BW/a +i+^ /2 (t))) 

X s (0) + 16u,\\v\\ LlB .+t + 16u;s%\\ LtB ,^ +e 2p \\Vf\\ L}B .-^ 

(3.21) +Ca\n(e- 1 )e*(x N/2 (t)\\v\\ LrB .-i + £ f '\\v\\ L]B ^\\v\\ L?B s^) ■ 
where 

(3.22) X s (0) = 4e||r || Bs + £ 2/3 (||r - 2wD[«o]||bj + ||r - 2wD[«o]|| B .-a) . 
Assuming that T* < To and noticing that 

(3.23) IMIri rN/2+i — ► as e — ► , 

we deduce from (|3.21[) - (|3.22[) and the continuity of t i— » X^/2(t) that there exists 
£o = £o(N, ||to|| b n/2,7-'/, uq) such that for any < e < Eq and any < t < T* , 

X N/2 (t) < A(e) 

with A(e) \ as e — > 0. This contradicts (|3.1[) of Theorem 13.11 and thus ensures that 
T* > T . |TTT6)) and (fl~18)) follow as well. To prove (fTTH]) we notice that from this 
last limit and ([OS]), \\t - 2wD[w]|| Ll b n, 2 -> 0, e 20 ' 1 ^ - 2wD[u]|| a b jv/2 -> and 

£ 2/3-i|| T _ 2wD[u]|| rA r n/2-2 -» 0. This gives the result since 2j3 - 1 + 2a < and 

thus 



I Li, B , 

r Je 



\t - 2uT>[ U ] || z i B . < £- 2Q ||r - 2wDM fl . a < e 2 ^ 1 ||r - 2 W D[ M ] \\ Ll B , 

J J e Tq J e Tq J e 
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Finally, for s > N/2, the proof follows the same lines using that 



(3.24) 




.N/2 + 1 



< e <»(— JV/2)|| w j| 




and thus with 



e = e (N, ||t || b jv/2, \\u \\ b n/2-i, \\v\\ 




This completes the proof of Theorems 11.21 and 11.31 
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